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1 Motivation 

Neutrino interactions in strong magnetic fields have gained a lot of attention be- 
cause of the problem of neutrino emission from a proto-neutron star. It has been 
argued that a small asymmetry in the neutrino emission can explain the high 
peculiar velocities of pulsars. Based on earlier workQ on neutrino dispersion 
relations in an external magnetic field, it has been shown[^ that neutrino oscil- 
lations can probably account for the pulsar kicks. On the other hand, there is 
a parallel stream of argument with the idea that neutrino opacities may change 
appreciably in magnetic fields. Asymmetric opacities can presumably account 
for the kicks. The opacity has recently been calculated by Roulet|3|, assuming 
the matrix element for the interactions remain unchanged in a magnetic field. 
He found that the opacities do change appreciably compared to the no-field 
case. However, he did not find any asymmetry vis-a-vis the direction of the 
magnetic field. The aim of the present work is to redo these calculations from 
first principles and show that the relevant cross-section is indeed asymmetric. 

2 Solutions of the Dirac equation in a uniform 
magnetic field 

We take the magnetic field B in the z-direction, and choose 



The Dirac equation in this field can be exactly solved. The energy eigenvalues 
are given by 



AQ=Ay=A^:^0, 



= -By. 



(1) 



El = +pI + 2neB , 



(2) 
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where n is a non-negative integer signifying the Landau level. All n > levels are 
doubly degenerate, whereas the n — level is non-degenerate. The eigenstates 
are of the form 



(3) 



where the notation p ■ Xi^^ stands for the dot product, setting the y-components 
equal to zero. Ug is the spinor, given by 



U+ = 







( \ 

Ma / (p\ 



(4) 



We have used the shorthand M„ = \/2neB, and a dimensionless variable 
The function appearing in Eq. is given by 



(5) 



(6) 



where Hn are Hermite polynomials, and Nn is a normalization which can be 
chosen arbitrarily. We choose 



Nr., 



'eB 



1/2 



(7) 



For the sake of consistency, we will define /_i = so that the solution 
vanishes for the zeroth Landau level. From the spinor solutions, the spin sum 
can be calculated: 

Pu{y,y*,n,i\) = ^ Us{y,n,p^)Us{yi,,n,p^) 

s 

+ {m(l-(T,)-t-;5|| -f^||75}/„_l(0/n-l(e*) 

+ M„{(7i-t-Z72)/„(0^n-l(e0 + (7l-*72)/n-l(0^n(^*)} , («) 

using the notations — j^.p'^^, — Jp.p'^^, where 

p|[ = (po,0,0,p,), = (p„0,0,po)- (9) 
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3 The fermion field operator 

The fermion field operator can be written as 



EE 



dpx dpz En + m 



2n 



2E„ 



s=± n=0 



(10) 



wliere /s(n,^^)is annihilation operator for the fermion, and /j^(n,p^) is the cre- 
ation operator for the antifermion in the n-th Landau level with given values of 
Px and Pz- The creation and annihilation operator satisfy the anticommutation 
relations 



fAn,P,)jl,{n',T() ^ = 5ss'5nn'S{px-p'x)5{Pz -p'J, 



(11) 



etc. The prefactor appearing in Eq. (^0|) is determined from the fact that the 
field operator should satisfy the anticommutation rule 



for X^ = X^. The one-fermion states are defined by 



27r 



(12) 



(13) 



the normalization constant chosen so that these states are orthonormal in a box 
with sides Lx, Ly and Lz- Then 



(14) 



where ipu denotes the term in Eq. (|10[) that contains the J7-spinors 



4 Inverse beta-decay 

We now apply the above formalism for the inverse beta-decay process 

z.e(fc)+n(P)^p(P') + e(p(,n')- (15) 

Assume eB ^ m^, so that the magnetic field effects on the proton and neutron 
spinors can be ignored. 

The interaction Lagrangian is 



V2Gf3 ■0(6)7'" -^-V'K) ■0(p)7M(5y +5A75)'0(n) 



(16) 
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where Gp = GpcosOc, gv = 1 and qa — 1.26. For the S'-matrix element, this 
gives 

% = A/2G;3yd"^(e(p(,n')|?(,)7''i^K)|;^e(fc) 

X ^p(F') |?(p)7m(3v + gA75)V'(„) I n{P)) . (17) 

The hadronic part of the matrix element is obtained with usual rules, which 
gives 

„i(P'-P)-X 



(18) 



For the leptonic part, we need magnetic spinors for the electron. Using Eq. (p^, 
we obtain this matrix element to be 



-ik-X+ip' -Xi^ I p 



C^(e)(y,n-',H')7^iwK)(^) 



(19) 



Putting these back into Eq. (|T^) and performing the integrations over all co- 
ordinates except y, we obtain 



Sf, = i27rfS^iP + k-P'-p') 
where 



Er,' + m 



2ujV 2PoV 2P^V2E„'L^L^ 



1/2 



, (20) 



= y/2GfJ^U^p){P')Y'{9V + gAl5)u(n){P) 

X / dye^'?"^ [17(e) (y,n',f^)7^iu(^^)(fc) 



using the shorthand 



(21) 



(22) 



Note that Eq. (gO|) does not have a J-function corresponding to the conservation 
of y-component of momentum, because the latter is not a good quantum number 
in this problem. The transition rate can now be written as 



1 



\Sj^/T = -{2^f5^{P + k-P'-p' 



16^-' - ■ - - - W^LuPoP^En' 

Putting in now the differential phase space for final particles, which is 

^^dv' —dv' ^ (fiP' 
2-K 2-K (2^)3 ^ 

where V — LxLyLz, we can write the differential cross section as 

1 

da 



(23) 



(24) 



647r 



1 o , , IS E,„i + m 

■^3(p + fc-p'-p'^ " ^ 



uoPaPoEn 



'^^^Up'^dp'.d'P'. (25) 



V 
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The square of the matrix element, averaged over initial neutron spin, is given 

by 



Here, H'^'^ is the hadronic part, which is 



The leptonic part contains magnetic spinors: 



(26) 



The integrations over y and j/* can be exactly performed Q, yielding 



27r 



eB {En' + m) 



(27) 



(28) 



(29) 



where 



eB 
+2Mngqh 



In- 



In-1 



qy_ 

JIb 



'eB 



(30) 



We calculate the cross section in the rest frame of the neutron. The axes are 
chosen such that incoming neutrino is in the x-z plane. Also, we assume that in 
the energy regime of interest, |P'| <C Wp so that the proton is non-relativistic. 
Then we obtain, for n' ^ 0, 



2 IGttG^ rrinmp 



eB En 



{.gl + 35i)t-'Ao + {gl - g\)k,A, 



(31) 



We now put it into Eq. (p5|). Integrations over and get rid of the S- 
functions. Integration over p'^ gives a factor LyeB. 



da = 



G} 5{Q + uj - En') 



47r 



LuEn' 



ig^ + 3g\)LuAo + {gl.~g\)k,A, dP'dp',, (32) 



where Q — m„ — nip. Integration over Py can now be performed using the 
orthogonality of the Hermite polynomials. This ensures that the M„-tcrm of 
Eq. ( po[) vanishes. The other two give 



d(Tn'>0 — 



eBGj 



27r 



-d{Q+LO-En') 



(g'v + SgD + ig'y-gl)^ 



dp', . (33) 
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The argument of the remaining (5-function vanishes for two different values of 
which are equal and opposite. The term proportional to gy — g\ then cancels 
out and we obtain 



O'n'>0 



35l 



Q + uj 



(34) 



For n' = 0, however, even the second term of Eq. ( |30| ) vanishes, so we need 
to redo everything beginning from Eq. (|3l|). The final result is: 



eBGl 
2tt 



fv + i9\)-{9l-9\)- 



(35) 
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Figure 1: Total cross section as a function of the magnetic field, normalized to 
the cross section in the field-free case. The initial neutrino energy is lOMeV. 
The solid and the dashed lines are for the initial neutrino momentum parallel 
and antiparallel to the magnetic field. 

For any given energy of the initial states, there is a maximum possible value 
for n', dictated by Eq. (||). The total cross section is obtained by performing 
the sum over all possible n' . We can then summarize our results as follows: 

• (To depends on the direction of neutrino momentum because of the 
term involving fc^. 

• Total a always contains (Toi s^nd is hence asymmetric. 

• The fractional amount of asymmetry depends on n'^^^. It is max- 
imum if n[^^^^ ~ 0, and decreases at smaller values of B as more 
Landau levels contribute. 

In Fig. |l|, we show the cross section as a function of B for neutrinos incident 
along and opposite to the direction of the magnetic field. For small values of 
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B, the distinction between the two directions is not seen at the scale of the 
plot. However, when B is so large that n'^^^^ = 0, the asymmetry is large. It is, 
in fact, as large as 18% between these two extreme directions. It seems highly 
plausible that such a large asymmetry in the cross section would produce enough 
momentum asymmetry to explain the pulsar kicks. These calculations were not 
finished at the time of the conference and therefore could not be presented [^. 

There is however another important issue to be settled before the momentum 
asymmetry is calculated. Neutrinos undergo many other reactions inside a dense 
star. The inverse beta decay has the largest cross section for the B = Q case 
of course, but there is no guarantee that other cross sections are not larger in 
strong magnetic fields. Our preliminary calculations for the neutrino-electron 
elastic scattering show that it can increase by several orders of magnitude in 
a magnetic field. If that is the case, then the opacities have to be calculated 
taking these processes into account. 
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